Pressure-induced structural, electronic, and magnetic effects in BiFeC>3 
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We present a first-principles study of multiferroic BiFeOa at high pressures. Our work reveals the 
main structural (change in Bi's coordination and loss of ferroelectricity), electronic (spin crossover 
and metallization), and magnetic (loss of order) effects favored by compression and how they are 
connected. Our results are consistent with the striking manifold transition observed experimentally 
by Gavriliuk et al. [Phys. Rev. B 77, 155112 (2008)] and provide an explanation for it. 
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Room-temperature multiferroic BiFeC>3 (BFO) is one 
of the most intensively studied materials of the moment. 
BFO is among the most promising multifcrroics from the 
applications perspective, and the latest results regarding 
its fundamental properties [lj and the engineering possi- 
bilities it offers [2] continue to fuel the interest in it. 

Indeed, recent works suggest the correlations between 
the structural, electronic, and magnetic properties of 
BiFeG"3 are not understood yet. We have evidence for 
a large sensitivity of BFO's conductivity to magnetic or- 
der [l], a metal- insulator (MI) transition driven by 
structural changes at high temperatures 0], metallicity 
at ferroelectric domain walls |5(, and even a spin-glass 
phase below 150 K [6J]. Yet, the results most revealing of 
the complex interactions in BFO may be those of Gavril- 
iuk et al. [7], |8( . These authors observed a pressure-driven 
diffuse transition, occurring at room temperature in the 
40-50 GPa range, that involves a structural change, loss 
of magnetic order, and metallization. The driving force 
behind such transformations, tentatively attributed to a 
spin crossover of Fe 3+ , remains to be clarified. 

Here we report on a first-principles study of the high- 
pressure behavior of BFO in the limit of very low tem- 
peratures (nominally, K). Our results reproduce the 
essential observations of Gavriliuk et al, thus reveal- 
ing the mechanisms that can lead to a manifold (struc- 
tural/electronic/magnetic) transition in this material. 

Methods- We used the local density approximation 
(LDA) [g] to density functional theory as implemented 
in the code VASP 10], including the so-called LDA+U 
correction of Dudarev et al. [ll| for a better treatment 
of iron's 3d electrons. The technicalities of our calcula- 
tions [12| are standard, but our use of the LDA+U asks 
for a comment: In this work we compared the energies of 
different electronic phases (insulating/metallic) in which 
the Fe 3+ ions display different spin states. Doing this 
accurately constitutes a challenge for any ab initio ap- 
proach; in particular, while commonly used to study such 
problems, the computationally-efficient LDA+U should 
be employed with caution in this context. For this rea- 
son, we repeated all our calculations for different values 
of the U parameter in the 0-4 eV range, and thus made 
sure our qualitative conclusions are reliable. (Unless oth- 



erwise indicated, the reported results are for U = 3 eV.) 

Phase transitions in BFO under moderately high pres- 
sures have been studied theoretically 0. In this work, 
however, we were not interested in the relatively mild 
effects so far investigated. Rather, we wanted to deter- 
mine whether pressure may induce profound changes in 
the electronic structure of the compound. Thus, we re- 
stricted our simulations to the 10-atom cell of the i?3c 
phase of BFO stable at ambient conditions (see inset in 
Fig- [J), assuming this is enough to capture the phenom- 
ena of interest. Then, as a function of volume (i.e., pres- 
sure), we performed structural relaxations for a variety 
of atomic, electronic, and magnetic configurations. That 
allowed us to identify a large number of possible phases 
and determine their relative stabilit y and p roperties. As 
done in other first-principles studies IH, 14 1 , we neglected 
the long-period spin cycloid in BFO. 

Structural and spin transitions.- Figure [T] displays our 
results for the equation of state of the BFO phases we 
found to be stable in some pressure range. The obtained 
pressure-driven transitions are better visualized in Fig. [2j 
which shows the pressure dependence of the key proper- 
ties. At ambient and moderately high pressures we ob- 
tained the usual ferroelectric i?3c phase of BFO, with 
G-type AFM order and Fe 3+ in a high-spin (HS) state. 
Then, we found that at about 36 GPa BFO undergoes a 
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FIG. 1: (color online). E(V) curves for the BFO phases here 
considered (see text). Inset: structure of the ferroelectric i?3c 
phase of BFO stable at ambient conditions. 
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FIG. 2: (color online). Pressure-dependence of the enthalpy 
(a), volume (b), structural distortions with respect to the 
ideal cubic perovskite structure (c), Fe's localized magnetic 
moment (d), and electronic band gap (e) for the phases of 
BFO here considered (see text). Dashed vertical lines mark 
the transition pressures. Bigger symbols indicate which phase 
is stable. Panel (a): We take as the zero of enthalpy the 
pressure-dependent result for the HS phase. Panel (c): For all 
phases, the upper (resp. lower) symbols correspond to the g- 
type (resp. u-type) distortions (see text); distortions in units 
relative to the pressure-dependent lattice vectors. Panel (e): 
The highlighted area is suggestive of the subtleties pertaining 
to the value of -E gap in the LS phase (see text). 



first-order phase transition to a phase with the Fe 3+ ions 
in a low-spin (LS) configuration. The concurrent drops 
in volume and iron's magnetic moment across the HS-LS 
transition can be seen, respectively, in Figs, [^b and[2JI. 
(The magnetic moments are nearly pressure-independent 
within both the HS and LS phases.) Finally, a metallic 
phase with no localized magnetic moments, denoted as 
"NS phase" where NS stands for "No Spin" , becomes sta- 
ble above 72 GPa. Such a paramagnetic metallic phase is 
typical of transition-metal oxides at high pressures 
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The main structural features of the obtained stable 
phases can be described in terms of g-type and u-type 
distortions of the ideal cubic perovskite structure, where 
g (resp. u) stands for gerade or even under inversion 
(resp. ungerade or odd under inversion). The evolu- 
tion with pressure of the distortions thus quantified is 
shown in Fig. For all the phases, the g-type distor- 
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HS/0GPa 2.31 2.42 12.9 
HS/36GPa 2.35 2.29 6.3 
LS/36GPa 2.51 2.25 0.0 



FIG. 3: (color online). Left: surroundings of Bi in the BFO 
structure. The six nearest-neighboring oxygens are split in 
two groups, denoted O(l) and 0(2) and colored differently, 
each of which is composed of three symmetry-related atoms 
that form a plane perpendicular to the [111] direction. Right: 
Structural parameters relevant to the HS-LS transition. 



tions are essentially oxygen-octahedron rotations as those 
occurring in the ferroelectric i?3c phase of BFO at ambi- 
ent pressure; our calculations show such rotations remain 
present under compression. The u-type displacements in 
the HS phase correspond to the usual ferroelectric distor- 
tion of BFO, dominated by the stereochemical activity of 
bismuth: the Bi atoms move along the [111] direction to 
approach the three O atoms forming a face of the neigh- 
boring oxygen octahedron. In the LS phase, though, we 
find a null u-type distortion (i?3c space group), which 
reflects a change in bismuth's coordination. As shown 
in Fig. [31 the LS phase presents Bi03 planar groups in 
which the three oxygens binding with one Bi atom be- 
long to three different 06 octahedra. Such planar groups 
are also present in the NS phase, but there they co-exist 
with a u-type displacement, along the rhombohedral axis 
(i?3c space group), in which the Bi and Fe atoms have a 
quantitatively similar participation. 

Thus, we found pressure favors two main modifications 
of the usual i?3c phase, namely, a change in Bi's coordi- 
nation and a low- or null-spin configuration of the Fe 3+ 
ions. We also checked such variations can exist indepen- 
dently (e.g., LS Fe 3+ can occur in absence of planar Bi03 
groups), which gives raise to various metastable phases. 
Indeed, we obtained many metastable phases (not shown 
here) with different atomic and electronic structures, in- 
cluding the occurrence of intermediate-spin Fe 3+ . We 
even found a mixed-spin phase, with one HS iron and 
one LS iron in our 10-atom cell, that becomes stable in 
the 35-37 GPa pressure range. 

Metallization at the HS-LS transition - As expected, 
for the HS phase we obtained an AFM insulating ground 
state. The energy difference between the AFM and 
ferromagnetic (FM) configurations varies from about 
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FIG. 4: (color online). Electronic density of states of the FM 
and AFM configurations of the LS phase of BFO at 50 GPa. 
The shaded area indicates the result for one Fe 3 + ion with its 
local magnetic moment pointing up. 



0.27 eV/f.u. at GPa to about 0.65 eV/f.u. at 36 GPa, 
which is compatible with the experimentally observed 
high transition temperature (Tn c1 ~ 643 K at GPa) 
and its increase with pressure As shown in Fig. [2^, 
the HS phase is an insulator throughout its stability 
range, and its energy gap decreases with pressure. 

According to our results, the properties of the LS phase 
are more complex. Throughout its stability range, the 
AFM and FM magnetic arrangements are nearly degen- 
erate, never differing by more than 0.03 eV/f.u. (More- 
over, the ground state shifts from AFM to FM as pressure 
increases.) This implies that the magnetic ordering tem- 
perature for the LS phase will be much (about 10 times) 
smaller than that of the HS phase. We also found the 
electronic structure of the LS phase depends strongly on 
the magnetic order. Figure [4] shows illustrative results at 
50 GPa: The AFM case presents a gap of about 0.8 eV, 
while for the FM order we get a half-metallic solution. 

These results imply that, if we were able to heat the 
LS phase up to room temperature (T room ), we would ob- 
tain a paramagnetic state. Further, this magnetically- 
disordered phase will appear as being metallic, since the 
the thermally-averaged equilibrium state will present a 
significant electronic density of states at the Fermi level. 
Thus, our calculations suggest that, at T room , BFO might 
undergo a pressure-driven HS-LS transition that brings 
about simultaneous structural, magnetic, and electronic 
(insulator-to-metal) transformations. 

Discussion - The picture of the HS-LS transition that 
emerges from our simulations is essentially identical to 
the one suggested by Gavriliuk et al. 0] to explain their 
observed pressure-driven transition. One should keep in 
mind, though, that the experiments of Ref. were per- 
formed at T room ; thus, the identification between the ex- 
perimental transition and our HS-LS transition relies on 
the assumption that our LS phase is the equilibrium state 



at -Troom ■ 

To resolve this issue, one would need to compute the 
temperature-pressure (T-P) phase diagram of BFO ab 
initio, which falls beyond the scope of this work. Yet, we 
can estimate the temperature stability range of the low-T 
phases by calculating their enthalpy difference with the 
cubic phase, as such a quantity should be roughly propor- 
tional to the temperature at which the transition to the 
cubic phase occurs. At GPa we obtained 0.88 eV/f.u. 
for the difference between the R3c HS phase and the 
lowest-lying cubic phase (which presents Fe 3+ in a HS 
state), a large value consistent with the experimen tally 
observed high transition temperature (1200 K) 0, Il6| - 
This enthalpy difference decreases gently with pressure: 
For example, at 50 GPa, the i?3c LS phase and the 
lowest-lying cubic phase (which presents Fe 3+ in a LS 
state) differ by about 0.60 eV/f.u. Hence, our calcula- 
tions suggest that, within the stability range of the LS 
phase, the transition into a LS cubic phase will occur 
at temperatures well above T room (roughly, in the 700- 
800 K range). Our results thus support the hypothesis of 
Gavriliuk et al. that the metallic phase observed experi- 
mentally at high pressures and T room contains LS Fe 3+ . 

We found additional support for this identification. 
The computed HS-LS transition pressure (about 36 GPa) 
agrees reasonably well with the observed one (40- 
50 GPa), and the wealth of competing (meta)stable 
phases that we found is consistent with the diffuseness 
of the experimental transition. Further, the results for 
V(P) in Fig. [2h are in good qualitative agreement with 
the experimental data, and the computed bulk modulus 
of the HS and LS phases are markedly different - the HS 
phase being significantly softer - as observed experimen- 
tally for the low- and high-pressure phases lTf . Finally, 



the pressure dependence of the enthalpy difference be- 
tween the low- and high-pressure phases was experimen- 
tally determined to be 12 meV/GPa at the transition 
region, and we obtain about 15 meV/GPa. 

It is not our purpose here to give a detailed elec- 
tronic picture of the HS-LS transition that we found. 
Let us just note our results for BFO strongly resem- 
ble what occurs in hematite (Fe20a), which undergoes 
a HS-LS transition with accompanying metallization at 
50 GPa [3 [H [U. Indeed, the conclusion of Ref. M 
that there is an enhancement of the metallic character of 
hematite's LS phase under pressure, driven by the broad- 
ening of the not-fully-occupied ti g bands of Fe 3+ , seems 
consistent with our findings. We should also note Gavril- 
iuk et al. have recently described the pressure-induced 
metallization in BFO in terms of a Mott-Hubbard pic- 
ture [1] . We have doubts about this interpretation: First- 
principles calculations show that BFO displays broad, 
significantly hybridized, valence bands (see Refs. O and 
l2ll for the usual HS phase of BFO and Fig. |4] for our 
LS phase), which suggests it may not be adequate to 
place this material on the Mott-Hubbard side of the usual 



4 



Zaanen-Sawatzky- Allen diagram [22j |. 

Our results allow us to discuss currently debated fea- 
tures of the T-P phase diagram of BFO. Following the 
discovery of the high-temperature metallic phase men- 
tioned above which has the ideal cubic perovskite 
structure with a 5-atom unit cell, it has been proposed 
this phase might extend its stabilityrange down to low 
temperatures and high pressures [J, (HI- Our results sug- 
gest to the contrary. While restricted to a 10-atom cell, 
our simulations do include the ideal cubic perovskite as a 
possible solution, and show this phase does not become 
the ground state under compression. For pressures ex- 
tending up to 100 GPa, we always observe symmetry- 
lowering distortions, associated to oxygen-octahedron 
tiltings and the stereochemical activity of Bi. This pre- 
diction is consistent with the experimental results of 
Gravriliuk et al. 0| and Haumont et al. 24 1 , who observe 
a non-cubic phase at high pressures and room tempera- 
ture. 

Finally, let us comment on the status of our quanti- 
tative results. As mentioned above, we repeated all our 
calculations for different values of the U parameter that 
defines the LDA+U functional. For the U values typi- 
cally used in studies of BFO 0, Q Hj| and other oxides 
with Fe 3+ [HI (i.e., in the 3-4 eV range), the obtained 
qualitative results are identical. At the quantitative level, 
the main difference is a positive shift of the transition 
pressures as U increases. For example, for U =4 eV the 
HS-LS transition occurs at 42 GPa, as compared with 
36 GPa for £7=3 eV. The differences are greater for the 
LS-NS transition, which occurs at about 130 GPa for 
[7=4 eV. This is not surprising, as a bigger value of U 
will tend to favor more insulating solutions. Hence, we 
think we can take our quantitative results as quite ap- 
proximate in what regards the HS-LS transition and the 
properties of the HS and LS phases. The LS-NS transi- 
tion pressure is, obviously, not well determined. 

In summary, we have revealed and explained the oc- 
currence of a manifold (structural/electronic/magnetic) 
phase transition in BiFeOa at pressures of about 40 GPa. 
We hope our results will contribute to a better under- 
standing of BiFeOs's T-P phase diagram and the com- 
plex interactions at work in this material. 
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